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We present a theory of quasiparticle Hall transport in strongly type-II superconductors within their 
vortex state. We establish the existence of integer quantum spin Hall effect in clean unconventional 
d^2_y2 superconductors in the vortex state from a general analysis of the Bogoliubov-de Gennes 
equation. The spin Hall conductivity a%y is shown to be quantized in units of This result 
does not rest on linearization of the BdG equations around Dirac nodes and therefore includes 
inter-nodal physics in its entirety. In addition, this result holds for a generic inversion-symmetric 
lattice of vortices as long as the magnetic field B satisfies Hd <^ B <^ Hc2- We then derive the 
Wiedemann-Franz law for the spin and thermal Hall conductivity in the vortex state. In the limit 
of T 0, the thermal Hall conductivity satisfies K^y = ^^{^)^Ta^y. The transitions between 
different quantized values of a^y as well as relation to conventional superconductors are discussed. 
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I. INTRODUCTION 

One of the fundamental characteristics of high tem- 
perature superconductors (HTS) is the apparent appli- 
cabihty of the d-wave BCS based phenomenology to 
the broad range of quasiparticle properties in the super- 
conducting state 1^. This is far from trivial property 
for materials known to exhibit strong electron correla- 
tions. Although the horizon is still not entirely clear 
and there remain few unresolved issues, examples being 
the temperature dependence of quasiparticle lifetimes or 
the penetration depth in the under-doped regime, [^|J|] 
the cumulative weight of evidence indicates that the low 
energy properties of cuprate superconductors are indeed 
governed by nodal quasiparticles with Dirac like disper- 
sion, as seen in assorted spectroscopic Q and transport 
measurements |^ . These and other experiments serve as 
the foundation for the "BCS-like d-wave paradigm" as 
it is currently used in both theory and interpretation of 
experiments. 

A natural question is how does this picture hold to- 
gether in a mixed phase, in presence of an external 
magnetic field and an array of superconducting vortices 
and, if it does, are there some special features of the d- 
wave quasiparticle phenomenology which could be used 
to deepen our understanding of high temperature su- 
perconductivity. Recent activity on the experimental 
front appears most encouraging in this regard. In par- 
ticular, measurements of the thermal Hall conductivity 
Kxy in cuprate superconductors are currently viewed as 
especially informative probes of quasiparticle dynamics. 
These measurements provide a clean way of extracting 
quasiparticle contribution to K^y as phonons, the other 
source of significant thermal conduction, do not couple 
to the magnetic field by virtue of being neutral. In addi- 
tion, contrary to what takes place in an ordinary electri- 



cal Hall conductivity experiment, vortices do not experi- 
ence strong Lorenz force since there is only heat current 
and no net electrical supercurrent to which vortices are 
strongly coupled. Thus, vortices tend to remain station- 
ary and their transport does not serve as a significant 
channel for heat conduction. 

Recently, measurements of K^y were conducted by 
Ong's group Q on YBCO samples with a very long mean 
free path. These experiments were carried out over a 
wide range of magnetic fields (up to ~ 14T) and at tem- 
peratures from T ~ 12.5K to above the superconducting 
transition Tc 90K. Unfortunately, the experiments 
are resolution limited below 12.5K as signal becomes too 
weak. At temperatures up to 25K and for a (T depen- 
dent) range of mag netic fields y/H/Tesla < T/25K the 
experiments seem to suggest a rather simple scaling form 
for Kr,y: 

Kxy{H, T) = const. X Vlf T. (1.1) 

On the theoretical front, the initial interest, largely in- 
spired by Gorkov and Schrieffer and, in a somewhat 
different context, by Anderson , was directed at the for- 
mation of "Dirac Landau levels" and their signatures in 
the quasiparticle thermodynamics and transport | p^ , pl| . 
This picture of Dirac Landau level quantization, while 
theoretically elegant and appealing, relies on the mini- 
mal coupling of the nodal BCS-like quasiparticles to the 
electromagnetic vector potential, A. The assumption of 
such minimal coupling seems innocent but it is not, on 
fundamental grounds. The physics behind the interac- 
tion of nodal quasiparticles with the external magnetic 
field and vortices was elucidated by Franz and Tesanovic 
]l2t . These authors devised a singular gauge transfor- 
mation which allows one to recast the original problem 
of BCS-like quasiparticles, moving under the combined 
influence of an external magnetic field and a superflow 
arising from vortex array, into that of Bloch particles 
moving in an effective non-uniform and periodic mag- 
netic field, the spatial average of which equals zero. This 
approach clearly demonstrates that the low energy por- 
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tion of the quasiparticle spectrum can be described as 
that of a relativistic Dirac particle minimally coupled to 
a fictitious U(l) gauge potential, i.e. the "Berry" gauge 
field a, which supplies the needed ±7r winding in the 
quantum mechanical phase of a quasiparticle as it encir- 
cles a vortex. Such half-flux Aharonov-Bohm scattering 
arises entirely through interaction of quasiparticles with 
vortices and it does not involve the external magnetic 
field explicitly. Thus, the cyclotron motion in a Dirac 
cone is caused exclusively by a time-reversal invariant 
"Berry" gauge field and cannot lead to any Dirac Lan- 
dau level quantization. Further progress came through 
the work of Marinelli, Halperin and Simon who ana- 
lyzed the quasiparticle excitation spectra of different vor- 
tex lattices and provided analytic symmetry arguments 
regarding the presence of nodal (zero energy) points in 
such spectra. These authors also devised a perturbation 
theory in the vicinity of nodal points which can be used 
to derive various results by analytic means. Various an- 
alytic results were also derived in the large anisotropy 
limit [^,^. Finally, the intricacies of Dirac equation 
in presence of half-flux Bohm-Aharonov scattering were 
addressed in Ref. ||f^, where the tight-binding regular- 
ization was introduced to restore the exact singular gauge 
symmetry in numerical calculations using the low energy 
(linearized) theory. Such regularization describes a lat- 
tice d-wave superconductor and is a natural choice for 
cuprates: after all, most of the microscopic theories of 
HTS start from tight-binding effective Hamiltonians. We 
should note that similar results for the quasiparticle spec- 
tra were also found in fully self-consistent calculations 
on the original Bogoliubov-deGennes (BdG) equations, 
using the basis formed by eigenstates of the magnetic 
translation group |f TUf^ ]. 

In this paper, we present a detailed analytical study of 
the quasiparticle Hall transport in a vortex state comple- 
mented by explicit numerical calculations |ig| ]. We con- 
sider a lattice d-wave superconductor of Ref. - This 
is important and necessary since the straightforward lin- 
earization of BdG equations drops curvature terms and 
results in n^y — Q,^. We employ the Franz and 
Tesanovic (FT) transformation so that we can use the 
familiar Bloch representation of the translation group in 
which the overall chirality of the problem vanishes. This 
should be contrasted with the original problem where 
the overall chirality is finite and the magnetic translation 
group states must be used instead. Naively, it might ap- 
pear that after an FT singular gauge transformation the 
effects of the magnetic field have somehow been trans- 
formed away since the new problem is found to have zero 
average effective magnetic field. Of course, this is not 
true. The presence of magnetic field in the original prob- 
lem reveals itself fully in the FT transformed quasipar- 
ticle wavef unctions. Alternatively, there is an "intrinsic" 
chirality imposed on the system which cannot be trans- 
formed away by a choice of the basis. One manifestation 
of this chirality is the Hall effect. The utility of the sin- 
gular gauge transformation in the calculation of the elec- 



trical Hall conductivity in the normal 2D electron gas 
in a (non-uniform) magnetic field was realized by Nielsen 
and Hedegard . They demonstrated that using singu- 
larly gauge transformed wavefunctions one still obtains 
the correct result, giving the electrical Hall conductance 
quantized in units of e'^/h if the chemical potential lies 
in the energy gap. In a superconductor, the question 
of Hall response becomes rather interesting as there is 
a strong mixing between particles and holes. Evidently, 
the electrical Hall response is very different from the nor- 
mal state, since charge is not conserved in the state with 
broken U(l) symmetry. Therefore, as pointed out in Ref. 
p2| , charge cannot be transported by diffusion. On the 
other hand, the spin is still a good quantum number p^ ] 
and it is natural to ask what is the spin Hall conductivity 
in the vortex state of an extreme type-II superconductor 
p3| . Moreover, every channel of spin conduction simulta- 
neously transports entropy [^,^^ and we would ex- 
pect some variation on Wiedemann-Franz law to hold 
between spin and thermal conductivity. 

As one of our main results, we derive the Wiedemann- 
Franz law connecting the spin and thermal Hall trans- 
port in the vortex state of a d-wave superconductor. In 
the process, we show that the spin Hall conductivity, a^y, 
just like the electrical Hall conductivity of a normal state 
in a magnetic field, is topological in nature and can be 
explicitly evaluated as a first Chern number characteriz- 
ing the eigenstates of our singularly gauge transformed 
problem [E5 23 2q]. Consequently, as T ^ 0, the spin 



Hall conductivity is quantized in the units of h/Sn when 
the energy spectrum is gapped, which, combined with 
the Wiedemann-Franz law, implies the quantization of 
Hxy/T. We then explicitly compute the quantized values 
of a^y for a sequence of gapped states using our lattice 
d-wave superconductor model in the case of an inversion- 
symmetric vortex lattice. Within this model one is nat- 
urally led to consider the BCS-Hofstadter problem: the 
BCS pairing problem defined on a uniformly frustrated 
tight-binding lattice. We find a sequence of plateau tran- 
sitions, separating gapped states characterized by differ- 
ent quantized values of a^y. At a plateau transition, level 
crossings take place and a^y changes by an even integer 
|p7| . Both the origin of the gaps in quasiparticle spectra 
and the sequence of values for a^y are rather different 
than in the normal state, i.e. in the standard Hofstadter 
problem |2^] . In a superconductor, the gaps are strongly 
affected by the pairing and the interactions of quasipar- 
ticles with a vortex array. The sequence of a'^y changes 
as a function of the pairing strength (and therefore in- 
teractions), measured by the maximum value of the gap 
function A |^ . Finally, we discuss the relation of our re- 
sults to those obtained within the continuum low energy 
(linearized) theory. 
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II. BOGOLIUBOV-DEGENNES HAMILTONIAN 



be written as a differential operator: 



The experimental evidence points towards well defined 
d-wave quasiparticles in cuprate superconductors in the 
absence of the external magnetic field. This suggests that 
to zeroth order fluctuations can be ignored and that one 
can think in terms of an effective BCS Hamiltonian, the 
simplest of which is written on the 2-D tight-binding lat- 
tice with the nearest neighbor interaction thus naturally 
implementing d^2_y2 pairing. In question is then the re- 
sponse of such a superconductor to an externally applied 
magnetic field B. All high temperature superconductors 
are extreme type-II forming a vortex state in a wide range 
of magnetic fields. This immediately sets up the contrast 
between B = and B 7^ situations: first, the problem 
is not spatially uniform and therefore momentum is not a 
good quantum number and second, the array of ^ vortex 
fluxes poses topological constraint on the quasi-particles 
encircling the vortices. Therefore, despite ignoring any 
fluctuations, the problem is far from trivial and demands 
careful examination. 

The natural starting point is therefore the mean-field 
BCS Hamiltonian written in second quantized form [pO[: 



i^ = /^xV4(x)(^(p-^Af 



/ij^a(x) + 

dx J dy[A(x,y)^|(x)i^|(y)-f A*(x,y)^4(y)VT(x)] 

(2.1) 



where A(x) is the vector potential associated with the 
uniform external magnetic field B, single electron energy 
is measured relative to the chemical potential /i, ■(/'q(x) is 
the fermion field operator with spin index a, and A(x, y) 
is the pairing field. For convenience we will define an 
integral operator A such that: 



AV'(x) 



dyA(x,y)V'(y). 



(2.2) 



In the strictest sense, on the mean field level this prob- 
lem must be solved self-consistently which renders any 
analytical solution virtually intractable. On the other 
hand, in the case at hand the vortex lattice is dilute for a 
wide range of magnetic fields, and by the very nature of 
cuprate superconductors having short coherence length, 
the size of the vortex core can be ignored relative to the 
distance between the vortices. Thus, to the first approxi- 
mation, all essential physics is captured by fixing the am- 
plitude of the order parameter A while allowing vortex 
defects in its phase. Moreover, on a tight-binding lattice 
the vortex flux is concentrated inside the plaquette and 
thus the length-scale associated with the core is implic- 
itly the lattice spacing S of the underlying tight-binding 
lattice. As shown in Ref. |16|, under these approxima- 
tions the d-wave pairing operator in the vortex state can 



A 



(x)/2 



(2.3) 



The sums are over nearest neighbors and on the square 
tight-binding lattice S — ±x, ±y; the vortex phase fields 
satisfy V x V0(x) = ^nzj^i <5(x ~ x^) with x^ denoting 
the vortex positions and (5(x— x^) being a 2D Dirac delta 
function; p is a momentum operator, and 



Vs = 



-1 



as 

if J 



±x 



(2.4) 



The operator rjs follows from the d-wave pairing: A = 
2AQ[cos{kxSx)—cos{kySy)]. For notational convenience we 
will use units where h — 1 and return to the conventional 
units when necessary. 

It is straightforward to derive the continuum version 
of the tight binding lattice operator A (see Ref. 
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^^,-R,{a.,A(x)}} 



{5„{5„A(x)}} + 

+-^A(x)[(a|</))-(a2^)], (2.5) 



but for convenience we will keep the lattice definition 
( |2.3D throughout. One can always define continuum as 
an appropriate limit of the tight-binding lattic e th eory. 
With the above definitions, the Hamiltonian ( ^.l[ ) can 
now be written in the Nambu formalism as 



H = J dx *^(x) Ho *(x) 



(2.6) 



where the Nambu spinor ^'^ — (jppipi) and the matrix 
differential operator 



h 

A* 



A 

-h* 



(2.7) 



In the continuum formulation h = Tr^(p — -A)^ — u, 
while on the tight-binding lattice: 



(2., 



t is the hopping constant and /i is the Fermi energy. The 
equations of motion of the Nambu fields ^ are then: 



ihi' = [^,H] = Ho^. 



(2.9) 



Note, that the Hamiltonian in Eq. (2.1) is our starting 
unperturbed Hamiltonian. In order to compute the lin- 
ear response to externally a ppli ed perturbations we will 
have to add terms to Eq. (2.1). In particular, we will 
consider two types of perturbations in the later sections: 
First, partly for theoretical convenience, we will consider 
a weak gradient of magnetic field (VB) on top of the 
uniform B already taken into account fully by Eq. (^.l|). 
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The VB term induces spin current in the superconductor 
pl| . The response is then characterized by spin conduc- 
tivity tensor which in general has non-zero off-diagonal 
components. Second, we consider perturbing the system 
by pseudo-gravitational field, which formally induces flow 
of energy (see [^2|-^) and allows us to compute thermal 
conductivity K^y via linear response. The advantage of 
these formal considerations are made clear in Section IV. 



A. Particle-Hole Symmetry 



magnetic unit cell 



The equations of motion ( |2.9| ) for stationary states lead 
to Bogoliubov-de Gennes equations Q 



(2.10) 
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The solution of these coupled differential equations are 
quasi-particle wavefunctions that are rank two spinors in 
the Nambu space, $^(r) = (u(r),w(r)). The single par- 
ticle excitations of the system are completely specified 
once the quasi-particle wavefunctions are given, and as 
discussed later, transverse transport coefficients can be 
computed solely on the basis of <l>'s. It is a general sym- 
metry of the BdG equations that if (u„(r), ti„(r)) is a 
solution with energy e„, then there is always another so- 
lution (— f * (r), u* (r)) with energy — e„ (see for example 
Ref. 13 ). 

In addition, on the tight-binding lattice, if th e ch emical 
potential /i = in the above BdG Hamiltonian (2/7), then 
there is a particle-hole symmetry in the following sense: 
if {un{v),Vn(y)) is a solution with energy e„, then there 
is always another solution e*'^''""="'"''*-'(u„(r), i;„(r)) with 
energy — e„. Thus we can choose: 



(2.11) 



where + (— ) corresponds to a solution with positive (neg- 
ative) energy eigenvalue. We will refer to this as particle 
hole transformation Ph- 



B. Franz-Tesanovic Transformation and Translation 
Symmetry 

In order to elucidate another important symmetry of 
the Hamiltonian {^^, we follow FT ||l2|,|6| and per- 
form a "bipartite" singular gauge transformation on the 
Bogoliubov-de Gennes Hamiltonian ( 2.10| ), 











(2.12) 



where (/>e(r) and (j)h{Y) are two auxiliary vortex phase 
functions satisfying 



(2.13) 



FIG. 1. Example of A and B sublattices for the square 
vortex arrangement. The underlying tight-binding lattice, on 
which the electrons and holes are allowed to move, is also 
indicated. 



This transformation eliminates the phase of the order pa- 
rameter from the pairing term of the Hamiltonian. The 
phase fields ^e(r) and ^/i(r) can be chosen in a way that 
avoids multiple valuedness of the wavefunctions. The way 
to accomplish this is to assign the singular part of the 
phase field generated by any given vortex to either (jje (r) 
or 0h(r), but not both. Physically, a vortex assigned to 
0e(r) will be seen by electrons and be invisible to holes, 
while vortex assigned to 0h(r) will be seen by holes and 
be invisible to electrons. For periodic Abrikosov vortex 
array, we implement the above transformation by divid- 
ing vortices into two groups A and B, positioned at {rf} 
and {rf } respectively (see Fig. |]). We then define two 
phase fields (j)^{r) and 4>^ [r) such that 

V X V(?;."(r) = 27rz^(5(r-rn, a = A,B, (2.14) 

i 

and identify 0e = 0'^ and cf)h — (f)^ ■ On the tight-binding 
lattice the transformed Hamiltonian becomes 



-|-criAo7?5e'-'^'- 



■P} (2.15) 



where 



1 e 

V = -V0 A; a : 

2 ^ c 



-(V(/.-4-V0^), 



(2.16) 



CTi and (T3 are Pauli matrices operating in Nambu space, 
and the sum is again over the n earest neighbors. Note 
that the integrand of Eq. (2.15) is proportional to the 
superfiuid velocities 



i^(V0"-fA) 

m* c 



a ■ 



A,B. 



(2.17) 
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and is therefore exphcitly gauge invariant as are the ofF- 
diagonal pairing terms. 

From the perspective of quasiparticles and vf can 
be thought of as effective vector potentials acting on 
electrons and holes respectively. Corresponding effec- 
tive magnetic field seen by the quasiparticles is B"fj = 
— ^ ^(V X v"), and can be expressed using Eqs. ( 2.14 ) 
and (2.15) as 



Ba 
off 



B 



t>oz 



(2.18) 



where B = V x A is the physical magnetic field and 
(f>o = hc/e is the flux quantum. We observe that quasi- 
electrons and quasi-holes propagate in the effective held 
which consists of (almost) uniform physical magnetic 
field B and an array of opposing delta function "spikes" 
of unit fluxes associated with vortex singularities. The 
latter are different for electrons and holes. As discussed 
in ||l^,0 this choice guarantees that the effective mag- 
netic fleld vanishes on average, i.e. (B"jj) = since we 
have precisely one ffux spike (of A and B type) per flux 
quantum of the physical magnetic fleld. Flux quanti- 
zation guarantees that the right hand side of Eq. (2. IS) 



vanishes when averaged over a vortex lattice unit cell con- 
taining two physical vortices. It also implies that there 
must be equal numbers of A and B vortices in the system. 

The essential advantage of the choice with vanishing 
(B^ff) is that and vf can be chosen periodic in space 
with periodicity of the magnetic unit cell containing an 
integer number of electronic flux quanta hc/e. Notice 
that vector potential of a field that does not vanish 
on average can never be periodic in space. Condition 
(B"jj) = is therefore cruc ial in this respect. The sin- 



gular gauge transformation (2.12) thus maps the original 
Hamiltonian of fermionic quasiparticles in finite magnetic 
field onto a new Hamiltonian which is formally in zero 
average field and has only "neutralized" singular phase 
windings in the off-diagonal components. 

The resulting new Hamiltonian now commutes with 
translations spanned by the magnetic unit cell i.e. 



[Tr, Hn] = 0, 



(2.19) 



where the translation operator Tr. = exp(iR-p). We can 
therefore label eigenstates with a "vortex" crystal mo- 
mentum quantum number k and use the familiar Bloch 
states as the natural basis for the eigen-problem. In par- 
ticular we seek the eigensolution of the BdG equation 
Hn'iP — £"0 in the Bloch form 

^„k(r) = e^''--$„k(r) = e"^ "" ^^^g 

where {Unk, Kik) are periodic on the corresponding unit 
cell, n is a band index and k is a crystal wave vector. 
Bloch wavefunction (r) satisfles the "off-diagonal" 
Bloch equation H{k.)^nk = ^nk^nk with the Hamilto- 
nian of the form 



(2.20) 



~zk r Tj ^«k r 



(2.21) 



Note, that the dependence on k, which is bounded to lie 
in the flrst Brillouin zone, is continuous. This will become 
important when topological properties of spin transport 
are discussed in the Section III C. 



III. SPIN CONDUCTIVITY 

Within the framework of linear-response theory 
spin dc conductivity can be related to the spin current- 
current retarded correlation function Dj^^^ through : 

(3.1) 

The retarded correlation function 1?^^, (il) can in turn be 
related to the the Matsubara flnite temperature correla- 
tion function 



/3 



e-"(r,j-(r)j-(0))dr 



(3.2) 



lim Dj^,{qi,q2,n)^D^,{in^n + tO). (3.3) 



In the Eq. (3^) the spatial average of the spin current 
j*(r) is implicit, since we are looking for dc response of 
spatially inhomogeneous system. In the next section we 
derive the spin current and evaluate the above formulae. 



A. Spin Current 

In order to flnd the dc spin conductivity, we must flrst 
flnd the spin current. More precisely, since we are looking 
only for the spatial average of the spin current j^(T) we 
just need its k — > component. In direct analogy with 
the B = situation |2^, we can deflne the spin current 
by the continuity equation: 



(3.4) 



where = §('0|j,V'T2: ~ V'jajV'ia;) is the spin density pro- 
jected onto z-axis. We can then use equations of motion 
for th e ^ fields (2.9) and compute the current density 
from (IJ). 

In the limit of g ^ the spin current can be written 
as (see Appendix B) 



(3.5) 



where the Nambu field = {i/j^^,tpi) and the general- 
ized velocity matrix operator satisfies the following 
commutator identity 



(3.6) 
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The equation ( |3.6D is a direct restatement of the fact 
that spin can be transported by diffusion, i.e. it is a 
good quantum number in a superconductor, and that 
the average velocity of its propagation is just the group 
velocity of the quantum mechanical wave. 

In the clean limit, the transverse spin conductivity a^y 
defined in Eq. (|3.l[) is 



Ai 



-trmrfxrnm 
^ y * X 



l0)2 



(3.7) 



where /,„ = (l + cxp(/3e„)) is the Fermi-Dirac distri- 
bution function evaluated at energy e,„. For details of 
the derivation see Appendix B. The indices m and n la- 
bel quantum numbers of particular states. The matrix 
elements F™" are 



(3.8) 



where the particle-hole wavef un ction s Um , Vn satisfy the 
Bogoliubov-deGennes equation (2.1C). Note that unlike 
the longitudinal dc conductivity, transverse conductivity 
is well defined even in the absence of impurity scatter- 
ing. This demonstrates the fact that the transverse con- 
ductivity is not dissipative in origin. Rather, as will be 
discussed in Section III C, its nature is topological. 



In the 
becomes 



limit of T — > the expression (3.7) for a 



Ai 



E 



T/mriT/nm 
*^ X * y 



*^ y ^ X 



(3.9) 



The summation extends over all states below and above 
the Fermi energy which, by the nature of the supercon- 
ductor, is automatically set to zero. 



Vanishing of the Spin Conductivity at Half 
Filling = 0) 



potential. Semiclassical analysis must then be started 
from this vantage point and the above conclusions are 
not straightforward, since the quasiparticle motion is ir- 
reducibly quantum mechanical. 

Here we present an argument for the full quantum me- 
chanical problem, without relying on the semiclassical 
analysis. We show that spin conductivity tensor (|3.7| ) 
vanishes at /x = due to particle hole symmetry ( p.ll[ ). 
First note that the Fermi-Dirac distribution function sat- 
isfies /(e) = 1 — /(— e). Therefore, the factor — /„ 
chang es sign under the particle-hole transformation Ph 
( ^.ll|) while the denominator (e^ — e„)^ clearly remains 
unchanged. In addition, each of the matrix elements V^"^ 

changes sign under Ph . Thus the double summation over 
all states in Eq. ( |3.7[ ) yields zero. 

Consequently the spin transport vanishes for a clean 
strongly type-II ECS d-wave superconductor on a tight 
binding lattice at half filling. Due to Wiedemann-Franz 
law, which we derive in the next section, thermal Hall 
conductivity also vanishes at half filling at sufhcicntly 
low temperatures. Note that this result is independent of 
the vortex arrangement i.e. it holds even for disordered 
vortex array and does not rely on any approximation 
regarding inter- or intra- nodal scattering. 



C. Topological Nature of Spin Hall Conductivity at 

T=0 

In order to elucidate the topological nature of a^y, we 
make use of the translational symmetry discussed in Sec- 
tion II B and formally assume that the vortex arrange- 
ment is periodic. However, the detailed nature of the 
vortex lattice will not be specified and thus any vortex ar- 
rangement is allowed within the magnetic unit cell. The 
conclusions we reach are therefore quite general. 

We will first rewrite the velocity matrix elements V^'" 
using the singularly gauge transformed basis as discussed 
in Section II B. Inserting unity in the form of the FT 
gauge transformation (2.12) 



It is useful to contrast the semiclassical approach with 
the full quantum mechanical treatment of transverse spin 
conductivity. In semiclassical analysis the starting un- 
perturbed Hamiltonian is usually defined in the absence 
of magnetic field B. One then assumes semiclassical dy- 
namics and no inter-band transitions. In this picture, if 
there is particle-hole symmetry in the original (B — 0) 
Hamiltonian, then there will be no transverse spin (ther- 
mal) transport, since the number of carriers with a given 
spin (energy) will be the same in opposite directions. In 
this context, similar argument was put forth in Ref. [Q. 
However, the problem of a d-wave superconductor is not 
so straightforward. As pointed out in Ref. in the 
nodal (Dirac fermion) approximation, the vector poten- 
tial is solely due to the superflow while the uniform mag- 
netic field enters as a Doppler shift i.e. Dirac scalar 



= {m\Vt,\n) = {m\U U-^V^^U U-^\n). (3.10) 



The transformed basis states U ^ | n) can now be written 
in the Bloch form as e^^''^\n-k) and therefore the matrix 
element becomes 



TT-ran 



= ( 



-ikr 



u 



(3.11) 



We used the same symbol k for both bra and ket be- 
cause the crystal momentum in the first Brillouin zone 
is conserved. The resulting velocity operator can now be 
simply expressed as 



1 9iJ(k) 



Ti dkn 



(3.12) 
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where -ffo(k) was defined in (2.21). Furthermore the ma- 
trix elements of the partial derivatives of H(k) can be 
simplified according to 

= -K-eZ%^\n,), (3.13) 



for m ^ n. Utilizing Eqs. (|1|) and (|l3|), Eq. (U) 



for crj,, can now be written as 



h 



e'"<0<£" 



dkv 



(^Fk)("k|7^I-)- (3-14) 



The identity J2e^<o<e^(\'^k)i'^k\ + hk)(".k|) = 1, can 
be further used to simplify the above expression to read 



'xy 



Stt 2TTi 



9mk 



dmi. 



dky 



9mk 



dky 



dm-L 



dkx 



(3.15) 



where crj^ is a contribution to the spin Hall conductance 
from a completely filled band m, well separated from the 
rest of the spectrum. Therefore the integral extends over 
the entire magnetic Brillouin zone that is topologically 
a two-torus T^. Let us define a vector field A in the 
magnetic Brillouin zone as 



i(k) = (mk|Vk|TOk), 



(3.16) 



where Vk is a gradient operator in the k space. From 
(3.15) this contribution becomes 



(3.17) 



where []z represents the third component of the vector. 
The topological aspects of the quantity in (3.17) were ex- 
tensively studied in the context of integer quantum Hall 
effect (see e.g. |3^) and it is a well known fact that 



1 



dk[Vk X i(k)], = Ci 



(3.18) 



where C'l is a first Chern number that is an integer. 
Therefore, a contribution of each filled band to a^y is 



-N 



(3.19) 



where N is an integer. The assumption that the band 
must be separated from the rest of the spectrum can 
be relaxed. If two or more fully filled bands cross each 
other the sum total of their contributions to spin Hall 
conductance is quantized even though nothing guaran- 
tees the quantization of the individual contributions. The 



quantization of the total spin Hall conductance requires 
a gap in the single particle spectrum at the Fermi energy. 
As discussed in the Section V, the general single parti- 
cle spectrum of the d-wave superconductor in the vortex 
state with inversion-symmetric vortex lattice is gapped 
and therefore the quantization of a'^y is guaranteed. 



IV. THERMAL CONDUCTIVITY 

Before discussing the nature of the quasi-particle spec- 
trum, we will establish a Wiedemann-Franz law between 
spin conductivity and thermal conductivity for a d-wave 
superconductor. This relation is naturally expected to 
hold for a very general system in which the quasi-particles 
form a degenerate assembly i.e. it holds even in the pres- 
ence of elastically scattering impurities. 

Following Luttinger |^^, and Smrcka and Stfeda 
we introduce a pseudo-gr avit ational potential 
X • g/c^ into the Hamiltonian (2.6) where g is a con- 



X 

stant vector. The purpose is to include a coupling to 

the energy density on the Hamiltonian level. This formal 
trick allows us to equate statistical (TV(1/T)) and me- 
chanical (g) forces so that the thermal current in the 
long wavelength limit given by 

j« = L«(r)(^rvi-Vx), (4.1) 

will vanish in equilibrium. Therefore it is enough to con- 
sider only the dynamical force g to calculate the phe- 
nomenological coefficient L'^^. Note that thermal con- 
ductivity Kxy is 



(4.2) 



When the BCS Hamiltonian H introduced in Eq.(2.1) 
becomes perturbed by the pseudo-gravitational field, the 
resulting Hamiltonian Ht has the form 



(4.3) 



where F incorporates the interaction with the perturbing 
field: 

F=^Jd^ «'t(x) (TJox + xHo) *(x). (4.4) 

Since x is a small perturbation, to the first order in x the 
Hamiltonian Ht can be written as 

HT = JdKil + |)*t(x) H, (1 -t- |)*(x) (4.5) 

i.e. the application of the pseudo-gravitational field re- 
sults in rescaling of the fermion operators: 



* ^ * = (1 |)^'. 



(4.6) 
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If we measure the energy relative to the Fermi level, the 
transport of heat is equivalent to the transport of energy. 
In analogy with the Section III A, we define the heat 
current j*^ through diffusion of the energy-density hx- 
From conservation of the energy-density the continuity 
equation follows 

V-j« = 0. 



In the limit of g ^ the thermal current is 



= - ( ^tv ^ _ vj/ V ^ 



(4.7) 



(4.8) 



For details see Appendix C. Note that the quantum sta- 
tistical average of the current has two contributions, both 
linear in x, 

(j?) = iiZ) + U?,)^ + )a.x. (4.9) 

The first term is the usual Kubo contribution to Lj^^ 
while the second term is related to magnetization of the 
sample |]38| for transverse components of n^y and van- 
ishes for the longitudinal components. In Appendix C 
we show that at T = the term related to magnetization 
cancels the Kubo term and therefore the transverse com- 
ponent of Kpjy is zero at T = 0. To obtain finite temper- 
ature response, we perform Sommerfeld expansion and 
derive Wiedemann-Franz law for spin and thermal Hall 
conductivity. 

As shown in the Appendix C 



where 



E 



X y 



jrmrt'W'nm 
y X 



(4.10) 



(4.11) 



Note that o'^y(^ = 0) = cr^^,{T = 0). For a superconduc- 
tor at low temperature the derivative of the Fermi-Dirac 
distribution function is 



dm 



6 



Substituting (4.12) into (4.10) we obtain 



(4.12) 



(4.13) 



where cr^^ is evaluated at T 
the hmit of T ^ 



0. Finally, using ( |4.2D , in 



47r2 
~3~ 



fee 



(4.14) 



We recognize the Wiedemann-Franz law for the spin and 
thermal conductivity in the above equation. As men- 
tioned, this relation is quite general in that it is indepen- 
dent of the spatial arrangement of the vortex array or 
elastic impurities. Thus, quantization of the transverse 
spin conductivity a^y implies quantization of Hxy/T in 
the hmit of T ^ 0. 



V. QUASIPARTICLE SPECTRUM AND 
QUANTIZED CONDUCTIVITY 

General features of the quasi-particle spectrum can 
be understood on the basis of symmetry alone. Since 
the time-reversal symmetry i s br oken, the Bogoliubov- 
de Gennes Hamiltonian Hq ( ^.1C| ) must be, in general, 
complex. According to the "non-crossing" theorem of 
von Neumann and Wigner js^, a complex Hamiltonian 
can have degenerate eigenvalues unrelated to symmetry 
only if there are three parameters which can be varied 
simultaneously. 

Since the system is two dimensional, with the vortices 
arranged on the lattice, there are two parameters in the 
Hamiltonian -ff(k) ( ^.21 ): vortex crystal momenta 



and ky which vary in the first Brillouin zone. There- 
fore, we should not expect any degeneracy to occur, in 
general^ unless there is some symmetry which protects 
it. Away from half-filling (/x ^ 0) and with unspecified 
arrangement of vortices in the magnetic unit cell there is 
not enough symmetry to cause degeneracy. There is only 
global Bogoliubov-de Gennes symmetry relating quasi- 
particle energy ek at some point k in the first Brillouin 
zone to — e_k. 

In order for every quasiparticle band to be either com- 
pletely below or completely above the Fermi energy, it 
is sufficient for the vortex lattice to have inversion sym- 
metry. This can be readily seen by the following argu- 
ment: Consider a vortex lattice with inversion symme- 
try. Then, by the very nature of the superconducting 
vortex carrying ^ fiux, there must be even number of 
vortices per magnetic unit cell and we are then free to 
choose Franz- Tesanovic labels A and B in such a way 
that v'^(— r) = — v^(r). To see this note that the ex- 
plicit form of the superfiuid velocities can be written as 
0: 



v"(r) = 



27rn 



^ik-(r-r°) 



(5.1) 



where a — A ox B and rf denotes the position of the 
vortex with label a. If the vortex lattice has inversion 
symmetry then for every vf there is a corresponding ^vf 
such that rf = — rf . Therefore, under space inversion I 



Iv^{r) = v^(-r) 



(5.2) 



Recall that the tight-binding lattice Bogoliubov- 
de Gennes Hamiltonian written in the Bloch basis ( 2.21 ) 
reads: 

s 

+aiAoVse'^^'^"^'''e''-^^+P^} (5.3) 

where 

v(r) ^ i(v^(r) + v^(r)); a(r) ^ i(v^(r) - v^(r)). 

(5.4) 
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As before cti and (73 are Pauli matrices operating in 
Nambu space and the sum is again over the nearest neigh- 
bors. It can be easily seen that upon applying the space 
inversion X to i?(k) followed by complex conjugation C 
and ia2 we have a symmetry that for every ek there is 
— ek, that is: 



-ia2CX -ff (k) XCia2 



(5.5) 



which holds for every point in the Brillouin zone. There- 
fore, in order for the spectrum not to be gapped, we 
would need band crossing at the Fermi level. But by 
the non-crossing theorem this cannot happen in general. 
Thus, the generic quasi-particle spectrum of inversion 
symmetric vortex lattice is gapped. As was established 
in the previous section, gapped quasi-particle spectrum 
implies quantization of the transverse spin conductivity 
a^y as well as Kxy/T for T sufSciently low. 

Precisely at half-filling (/i = 0) a^y must vanish on 
the basis of particle- hole symmetry (see Section III B). 
We can then vary the chemical potential so that 7^ 
and break particle-hole symmetry. Hence, the chemical 
potential /i can serve as the third parameter necessary for 
creating the accidental degeneracy, i.e. at some special 
values of pi* the gap at the Fermi level will close (see 
Fig. H). This results in a possibility of changing the 
quantized value of a^y by an integer in units ofh/Sn (Fig. 
g). [0 By the very nature of the superconducting state, 
we achieve plateau dependence on the chemical potential. 
This is to be contrasted to the plateaus in the "ordinary" 
integer quantum Hall effect which are due to the presence 
of disorder. In our case, the system is clean and the 
plateaus are due to the magnetic field induced gap and 
superconducting pairing. 

Similarly, we can change the strength of the electron- 
electron attraction, which is proportional to the maxi- 
mum value of the superconducting order parameter Aq 
while keeping the chemical potential fi fixed. Again, as 
can be seen in Fig. ^, at some special values Ag the 
spectrum is gapless and the quantized Hall conductance 
undergoes a transition. 



VI. SCALING FUNCTIONS 

As shown in the previous section, the quasi-particle 
spectrum of a d-wave superconductor in vortex state is 
gapped, in general. Therefore, at temperatures T which 
are much less than the A^ (magnetic field B induced 
gap), T-dependence of thermal conductivity be 
determined uniquely. Moreover, B-dependence of n^y 
comes entirely from the spin conductivity ctJj^(B). That 



47r2 
~3~ 



(6.1) 




0.4 0.6 
chemical potential n[t] 

FIG. 2. The mechanism for changing the quantized spin 
Hall conductivity is through exchanging the topological 
quanta via ( "accidental" ) gap closing. The upper panel dis- 
plays spin Hall conductivity a%y as a function of the chemical 
potential pL. The lower panel shows the magnetic field induced 
gap Am in the quasiparticle spectrum. Note that the change 
in the spin Hall conductivity occurs precisely at those val- 
ues of chemical potential at which the gap closes. Hence the 
mechanism behind the changes of a^y is the exchange of the 
topological quanta at the band crossings. The parameters for 
the above calculation were: square vortex lattice, magnetic 
length I = 45, A = O.lt or equivalently the Dirac anisotropy 
an = 10. 




FIG. 3. The upper panel displays spin Hall conductivity 
a^y as a function of the maximum superconducting order pa- 
rameter Aq. The lower panel shows the magnetic field induced 
gap in the quasiparticle spectrum. The change in the spin Hall 
conductivity occurs at those values of Aq at which the gap 
closes. The parameters for the above calculation were: square 
vortex lattice, magnetic length I = 45, = 2.2t. 
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Curiously, if the above equation is naively combined with 
Simon and Lee scaling 



(6.2) 



K.,y(T,B)=T^F,y[ V 



B 
f 



the scaling function would be determined up to a propor- 
tionality constant C: 



K,y{T,B)^CTVB. 



(6.3) 



In the recent experiments of Ong et. al. this is pre- 
cisely the scaling seen in the temperature range up to 
~ 25K. Although tempting, the above arguments are 
unlikely explanation for the experiments since the mea- 
surements are done at rather high temperatures. 



VII. CONTINUUM VS. LATTICE THEORY 

The previous discussion concentrated on the tight- 
binding formulation of the problem which is important 
if the magnetic field is relatively large and if there is a 
strong interaction between the underlying ionic lattice 
and the quasiparticles. In usual experimental situations, 
however, the magnetic length is much longer than the 
inter-ionic spacing and we would expect that the length- 
scale associated with the ionic lattice becomes unimpor- 
tant at low energies. This leads to a continuum formula- 
tion of the theory (see Section I). 



A. Linearized Continuum Theory 

On the basis of the B — problem, we expect that 
the low energy physics is confined in fc-space around four 
nodal points on the Fermi surface. In order to treat the 
effect of the magnetic field and of the vortex lattice on the 
low energy properties of the spectrum, we would ideally 
like to construct a Hamiltonian which treats each of the 
four nodes independently. A natural way to do this was 
suggested by Simon and Lee |Q and later extensively used 
by others [[2pjll,|l§. 

As was pointed out in Ref. there is a subtle 
problem associated with numerical implementation of 
the continuum linearized version of the theory. On the 
physical grounds any given problem involving charge e 
(quasi)particle interacting with a +^ vortex must lead 
to an identical eigenvalue spectrum as for a — ^ vortex as 
long as the wavefunction of the particle vanishes precisely 
at the vortex. This is an exact singular gauge symmetry 
of our physical problem and therefore it should be present 
at all stages of any approximation. Extensive numerical 
study of the linearized approximation fl^ pointed out 
that this symmetry is weakly violated in the quasipar- 
ticle spectra. The problem persists irrespective of the 
techniques chosen to diagonalize the linearized Hamilto- 
nian i.e. it is present in real space as well as momentum 



space representation. The issue was resolved in Ref. |T^] 
by regularizing the problem on the tight-binding lattice 
which offered many advantages, some of which are uti- 
lized in this paper. On the other hand, the concept of 
single node physics is lost in a tight binding formula- 
tion as the latter naturally describes the physics of all 
four nodes. It would be desirable to formulate a properly 
regularized Hamiltonian describing single node physics 
only. From the work of Franz and Tesanovic ]l^ , the 
most natural candidate is a free, anisotropic, massless 
Dirac Hamiltonian 



VfPx VAPy 
VAPy -Vppx 



(7.1) 



pertu rbed by scalar and vector potentials defined in Eq. 
(HI): 



-w' - ^ f ^^^^ 



syJ 



(7.2) 



vf and va are Fermi and gap velocities respectively [[l2| . 
We wish to argue that the above Hamiltonian is well 
defined and does not suffer from singular gauge asym- 
metry, provided that suitable boundary conditions and 
their transformations are also specified. The discrepan- 
cies encountered in the numerics are believed to result 
from neglecting the boundary conditions and the effect 
of the singular gauge transformation on them. 

In the context of high energy physics, a similar problem 
was studied extensively. It was found that the Dirac- 
type equations in the presence of a single Aharonov- 
Bohm string require self-adjoint extensions which enlarge 
the Hilbert space by the wavefunctions with r~2 radial 
behavior. In addition, the wavefunctions depend on a 
dimensionless parameter O which specifies the bound- 
ary conditions at the core. A consistent procedure en- 
abling one to construct the needed self-adjoint extension 
involves the theory of von Neumann deficiency indices 

The above considerations can be illustrated on a well 
studied problem of a Dirac particle in the field of a single 
string . We will restrict ourselves to the case of a 

massless Dirac particle and || flux carried by the string. 
The Hamiltonian 



(i ^~ef)ij} = 



(7.3) 



in this case is axially symmetric and after separation of 
variables 



V'(r, (j)) 



(7.4) 



one obtains the equation for the radial wavefunctions: 

i^.x.(_^(^«_,)-^(^;^))x(0^i.x(.) 



(7.5) 
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where = ^ + n is a half-integer. The solutions of this 
equation can be expressed through Bessel functions as 



Xr OC 



Je.i\E\r) 
iJei^+i)i\E\r) 



(7.6) 



where e = ±1. The square integrability condition deter- 
mines sign of e for all v except v = —1/2. Both choices 
of e result in a wavefunction diverging as ^1 \pr while still 
remaining square integrable. The requirement of reg- 
ularity for all the solutions at r = turns out to be 
too restrictive and results in numerous pathologies such 
as incompleteness of the basis. The problem is solved 
by using von-Neumann theory of self-adjoint extensions 
Since the radial B.^ operator in (7.f;) defined 
in the domain of regular functions has deficiency indices 
(1,1), one is forced to extend the Hilbert space by re- 
laxing the condition of regularity and allowing for the 
wavefunctions [H2| 



1 



i sin O 



cos 9 



(7.7) 



The angle 8 G (7i'/4, 7r/4 -I- tt), determined by the details 
of the vortex core, parameterizes the self-adjoint exten- 
sion. 

The energy eigenstates for v = —1/2 are given by 



sin^J-i/ad^k) + (-l)"cosMJi/2(|Sk) 
sinMJi/2(|^k) - (-l)"C0SMJ-l/2(|Sk) 



(7.8) 



where the parameter /i = if n is even, while for n odd 
/i = TT — 6. In addition, for some values of 8 there is a 
bound state [Q. 

We encounter a peculiar situati on i n which the long 
distance physics governed by Eq. ( |7.3| ) still depends on 
boundary conditions imposed at the core. Linearization 
of the original BdG problem neglects terms "small" in 
conventional sense which, neverthel ess, remain effectively 
present in 8 dependence of Eq. (7.8) . Although the 



linearized equation ( |7.3| ) has no length scale, the core 
physics, now effectively shrunk to a circle with vanish- 
ing radius around the singularity, manifests itself in the 
wavefunctions (7.7) with the long power-law tails. In 



general, there are only two values of corresponding to 
the pure 5-function magnetic field |4^ . However, in the 
case of other contact interactions present in the core, the 
parameter 8 can be arbitrary [[45| . 

Various regularizations of the problem were studied 
leading to different choices of the self-adjoint extension 
p6t . One possibility is to put an impenetrable cylinder 
of small but finite radius p around a vortex thus forcing 
the wavefunctions to vanish on the surface of the cylin- 
der. Such a boundary condition specification will im- 
mediately restore the ± || vortex invariance because the 
value of the wavefunction at the surface of the cylinder 
will not transform under singular gauge transformation. 
Of course, the true boundary conditions must descend 



from the original physical problem, in this case from the 
self-consistent solution of the full Bogoliubov-de Gennes 
equation including the core region. The quasi-particle 
wavefunctions do not, in general, vanish around the vor- 
tex pTj . Therefore, not only the quasi-particle wavefunc- 
tions but also the constraints imposed on them must be 
properly transformed under singular gauge transforma- 
tion. 



B. Beyond the Linearized Theory 

As pointed out by Simon and Lee and then later 
by Ye pO| ], linearized theory predicts vanishing trans- 
verse conductivities. In order to include the description 
of Hall effects, the curvature terms must be included. 
Without regularization, the perturbative analysis of the 
curvature terms are far from straightforward. This has to 
do with the fast rise of both wavefunctions and perturba- 
tions around the core. If the regularization is taken into 
account and the wavefunctions are forced to vanish at 
some radius p around the core, the new length scale p will 
break the scale invariance of the Dirac equation. More- 
over, the perturbative analysis of the continuum equation 
would predict that the contribution of each node to al.y 
can be either or ±i in units of fij^-K ||l^. Thus, if all 
four nodes are included, the only possible values of cr^y 
are or ±2. On the other hand, explicit evaluation of cr^j^ 
via the tight-binding lattice regularization indicates that 
a wide range of values for a^y is possible (see Figs. ^ and 
An interesting question remaining for future study 
is how to implement the picture presented in Ref. [ p^ , 
of individual Dirac nodes changing a%y by ±1 through 
action of the perturbatively determined mass term as an 
"elementary building block" of more complicated band 
crossings. 

One way to reach the continuum limit is to make the 
magnetic length / much longer than the tight binding 
lattice constant (5. The drawback associated with the 
tight-binding regularization of the problem stems from 
the fact that the continuum limit is not smooth and un- 
less an exact analytical solution is found it is complicated 
to analyze numerically. The similar situation occurs in a 
tight-binding Hofstadter problem of a 2D normal electron 
gas on the lattice. 



VIII. CONVENTIONAL SUPERCONDUCTORS 

Although the previous analysis was focused on uncon- 
ventional d-wave superconductors, the main results can 
be directly generalized to conventional 2D s-wave su- 
perconductors in the vortex state: the transverse spin 
Hall conductivity a%y is quantized in units of ?i/87r 

and by Wiedemann- Franz law n^y = "^"^^y 
as T — > 0. At low magnetic fields, the quasiparticle 
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spectrum is gapped by the virtue of the lowest Caroli- 
Matricon-de Gennes vortex core bound state having en- 
ergy ^ jtp. In the vortex lattice the CMdG bound 
states are extended and form bands. If the vortices are 
disordered, the states in the band tails will be localized 
and the effective quasiparticle gap will be increased be- 
cause the localized states do not contribute to transport. 
As the magnetic field is increased, the spectrum becomes 
progressively altered and at some critical field B* it de- 
velops nodal points (see Ref. @]). The transition be- 
tween states with two different quantized values of a%y 
happens at the B induced gap closings. Therefore, a se- 
ries of transitions between different spin Hall states is 
predicted to occur in the conventional superconductors 
as well. 



more complex. 

By varying some external parameter, for instance the 
strength of the electron-electron attraction which is pro- 
portional to Aq or the chemical potential /i, the gap 
closing is achieved i.e. = 0. The transition be- 

tween two different values of crj^^ occurs precisely when 
the gap closes and topological quanta are exchanged. 
The remarkable new feature is the plateau dependence 
of a^y on the Aq or It is qualitatively different from 
the plateaus in the ordinary integer quantum Hall effect 
which are essentially due to disorder. In superconduc- 
tors, the plateaus happen in a clean system because the 
gap in the quasiparticle spectrum is generated by the su- 
perconducting pairing interactions. 



IX. CONCLUSIONS 

In conclusion, we examined a general problem of 2D 
type-II superconductors in the vortex state with inver- 
sion symmetric vortex lattice. The single particle exci- 
tation spectrum is typically gapped and results in quan- 
tization of transverse spin conductivity a^y in units of 
h/Sn The topological nature of this phenomenon is 
discussed. The size of the magnetic field induced gap Am 
in unconventional d-wave superconductors is not univer- 
sal and in principle can be as large as several percent of 
the maximum superconducting gap Aq . By virtue of the 
Wiedemann-Franz law, which we derive for the d-wave 
Bogoliubov-de Gennes equation in the vortex state, the 

thermal conductivity K^y — ^ (^^) ^^^y as T ^ 0. 

Thus at T <C Am the quantization of Hxy/T will be ob- 
servable in clean samples with negligible Lande g fac- 
tor and with well ordered Abrikosov vortex lattice. In 
conventional superconductors, the size of Am is given by 
Matricon-Caroli-deGennes vortex bound states ~ A'^/ep- 
In real experimental situations, Lande g factor is not 
necessarily small. In fact it is close to 2 in cuprates . 
The Zeeman effect must therefore be included in the anal- 
ysis. Nevertheless, detailed numerical examination of the 
quasiparticle spectra reveals that the spectrum remains 
gapped for a wide range of physically realizable parame- 
ters even if g = 2. Thus, although the Zeeman splitting 
is a competing effect, in general it is not strong enough 
to prevent the quantization of a^y and consequently of 

^xy 1^ • 

We have explicitly evaluated the quantized values of 
al.y on the tight-binding lattice model of d^i-yi-viaYe 
BCS superconductor in the vortex state and showed that 
in principle a wide range of integer values can be ob- 
tained. This should be contrasted with the notion that 
the effect of a magnetic field on a d-wave superconductor 
is solely to generate a d + id state for the order param- 
eter, as in that case a^y = ±2 in units of h/Sn. In the 
presence of a vortex lattice, the situation appears to be 
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APPENDIX A: GREEN'S FUNCTIONS: 
DEFINITIONS AND IDENTITIES 

The one-particle Green's function matrix |Q is defined 

Ga/3(riTi;r2r2) = -(T,*„(l)^)5(2)) (Al) 

where Tr denotes imaginary time ordering opera- 
tor and a = 1,2 denotes components of a Nambu 
spmor ^'t(l) which is a shorthand for \E'^(ri,Ti) = 
(V'|(ri,Ti), Vi(ri,ri)) 



Due to the time ind epen dence of the Hamiltonian (2.1), 
the Green's function ([Al[ ) depends only on the imaginary 
time difference t = ti — r2. Therefore, its Fourier trans- 
form is given by 



G(ri,r2;zw) = /^^ e^"^G(ri, r2, T)dT 
G(ri,r2;T) = i ^ e"'"'G(ri, r2; zw). 



(A2) 



Here /? = l/{kBT), T is temperature, and the fermionic 
frequency uj ^ {21 + l)7r//3, I e Z. 

Using the above relations, it is straightforward to de- 
rive the spectral representation of the following correla- 
tion functions between ^' and its imaginary time deriva- 
tives 5,-* = 

f(T.vl,„(l)*^(2)) 

(r.^'a(l)*^(2)) ^1^ 

(T.vl/„(l)vi/t(2)) p2^^ 
[(X,vi,„(l)vi;t(2)) 



de 



A(ri,r2;e) 



-1 
-e 



(A3) 

The spectral function A(ri,r2;e) can be written in 
terms of eigenfunctions 4'n(r) — (u„(r), i;„(r))-^ of the 
Bogoliubov-deGennes Hamiltonian i^o in the form: 

i„^(ri,r2;e) :=^5(e-e„)$„„(ri)$J^^(r2), (A4) 



where e„ is the eigen-energy associated with an eigenstate 
labeled by the quantum number n. Substituting (A4) 
into (A3) we can write the above correlation functions 
solely in terms of the eigenfunctions $„ (r) : 



(T.vi,„(l)^-t(2)) 
(T.vl/„(l)vi/t (2)) 
[(X.vi;„(l)vi;t(2)) 



.*na(ri)$l^(r2) 



where VI = 27rfc/ /?, k € Z, is an outside bosonic fre- 
quency. The sum (A6) can be evaluated using standard 
techniques (see e.g. [|36t) and yields: 



fn fn 



ifl' 



(A7) 



Here /„ is a short hand for the Fcrmi-Dirac distribution 
function /(e„) = (l + exp(/3e„)) . 

In order to derive spin and thermal currents we will 
need the explicit form of the generalized velocity operator 
V introduced in (3.6): 



(A8) 



where the gap velocity operator is given by 
and the canonical momentum equals 



^ = _l5e^^^'+'(P-^^)"^' +h.c. 



(AlO) 



The following identities for operators va and tt will be 
used in the next section: 



7r**t.7r^' = -iW-(7r**t ^) + (7r*)2^rt 5» 

(All) 



^-tA^-A*^ * = -V-(*^ VA*-VA*^ *) (A12) 

The above equations are straightforward to derive in con- 
tinuum, while on the tight-binding lattice Eqs. (All) and 
(A12) imply a symmetric definition of the lattice diver- 
gence operator. 

The identities ( [All , Al^ ) explicitly ensure that the 
generalized velocity operator is Hermitian i.e it satis- 
fies the following identity: 



JdrV*¥{r) *(r) = j dvYl^-^liv) *„(r) 

= j dr^'t(r) V*(r). (A13) 



APPENDIX B: SPIN CURRENT AND SPIN 
CONDUCTIVITY TENSOR 



In the calculations that follow we will also encounter 
Matsubara summations over the fermionic frequencies 
w = (2/ + l)7r//3, ; e Z, of the form: 



-y 



1 



{iuj — e„)(ir2 + iuj — em) 



(A6) 



The time derivative of spin density ps = |(V'|V'T 
i/'I'f/'x) can be written in Nambu formalism as 



(Bl) 
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Using equations of motion (2.9) together with the exphcit 



form of Ho operator (2.7) we obtain 

^2 



2\ 2m 



2 / * \ 2 

2m 2m 



(B2) 



Using the identities (A11.A12) it is easy to show that 



(B3) 



wher e th e generalized velocity operator Vf^ is defined in 
Eq. ( A8 ) , and the last equality follows from the conti- 
nuity equation (3.4) relating the spin density p'^ and the 
spin current j**. Upon spatial averaging and utilizing Eq. 
( A13| ) we find the g ^ limit of the spin current 



(B4) 



APPENDIX C: THERMAL CURRENT AND 
THERMAL CONDUCTIVITY TENSOR 

In order to calculate thermal currents and thermal con- 
ductivity in the magnetic field we introduce a pseudo- 
gravitational potential x = r • g/c^ |]3^-^. This formal 
procedure is useful because it illustrates that the trans- 
verse thermal response is not given just by Kubo formula, 
but in addition it includes corrections related to magne- 
tization. Throughout this section h = 1. The pseudo- 
gravitational potential enters the Hamiltonian, up to lin- 
ear order in %, as 



Ht^ Jd^ [1 + |)vl/t(x) Ho (1 + |)*(x) 



(CI) 



where Ho is the Bogoliubov-deGennes Hamiltonian (^.7|). 
The equations of motion for the fields 5* thus become 



ii'^[^,HT] = (l + |)^o(l + |)* 



(l+x)^o4'-iV^xU^*. (C2) 



The ev alua tion of the spin current-current correlation 
function (3.2) is straightforward and yields: 



(r,*t(i)$(2)^t(3)$(4)) , Ut (B5) 

3, 4 ^ (y, 0) 

Here ^'(1) = ^'(ri,Ti) and similarly the operator U^(2) 
acts only on functions of V2. Using Wick's theorem, iden- 
tity (A5), and upon spatial averaging over x and y we 
obtain 



D^,{in) = — Y,{AVAm){m\VMSnmm- (B6) 



The last equality follows from the commutation relation 
K^. (Note that for x 7^ the Eq. (^) differs from Eq. 
(2.9). Throughout this section ^ will refer to Eq. (|C^ ) 
unless explicitly stated otherwise). 

To find thermal current j'^ we start with the continuity 
equation 



V • j'3 = 



(C3) 



The Hamiltonian density hr follows from Eq. (CI) and 
reads 

+ \ [^i/^cp^p + A^^^-l (C4) 



The double summation extends over the eigenstates \m) 
and \n) of t he Hamiltonian ( p.7[ ), and S',„„(if2) is given 
in Eq. (AT). Analytically continuing ifi — > il -I- ?0 we 
finally obtain the expression for the retarded correlation 
function: 



T/nm T/mn 



4 ^ €,r - e,- 



n + io 



ifn-frn). (B7) 



where is defined in (B^) and /„ is a short hand 
for the Fermi-Dirac distribution function /(£«) = (l + 

exp (/ ?£„) ) ^. Finally, we substitute the last equation 
into (3J^) to obtain 



+-2 T/nmT^mn 
^r..- = ^ E - frn) (B8) 



Ai ^ (e„ — 

ran ^ " 



where 5* = (1 -I- ^)^'. Taking the time deriv ative of the 
Hamiltonian density Kt and using Eqs.(All) and equa- 
tions of motion (C2) we obtain 



1 



it 



- ( *aAa;3^'/3 - *„Aa/3*/3 



where we introduced matrix operators 

Aq/3 = 



A 

A* 





< 



(C5) 



(C6) 



Here ^ = (1 + f )* and tt^ is defined in Eq. (|A1C|) 



Finally we use (A12) to extract j'^ from (C5). Upon 
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spatial averaging and using the Hermiticity of ( A13| ) 
the thermal current j'^ reads 



=U*^^M*-*V^* 



(C7) 



Note that the expression for j'^ contains two terms 



(C8) 



where Jq is independent of x ^-nd is linear in x- Ex- 
plicitly: 



(C9) 



and 



(CIO) 

where {a, 6} = a6 + 6a. Analogously to the situation in 
the normal metal, the thermal average of does not in 
general vanish in the presence of the magnetic field . 

The linear response of the system to the external per- 
turbation can be described by 

(j?) - + = -{K^u+M^.)d.x ^ -L%d.X 

(Cll) 



where 



S d^x 



= — lim 



(C12) 



is the standard Kubo formula for a dc response [ p6[ , 
P^(r2) being the retarded current-current correlation 
function, and 



Sd.X 

n n 

(C13) 

is a contribution from "diathermal" currents |^^. Note 
that the latter vanishes for the longitudinal response 
while it remains finite for the transverse response. As 
we will show later in this sectio n, at T = ther e is an 
important cancellation between (|C12| ) and (|C13| ) which 
renders the thermal conductivity k^i, well-behaved and 
prevents the singularity from the temperature denomina- 
tor in Eq. ( p^ . 

The retarded thermal current-current correlation func- 
tion P^(r2) can be expressed in terms of the Matsubara 
finite temperature correlation function 



Jo 



Trj^{r,r)Mr',0)^ (CM) 



as 



(C15) 



The current j'^(r, r) in Eq. (C14) is given by 
Mr) - \i^HT)V^dr^{T) - dr^Hr) V^^ (t)) (C16) 

As pointed out by Ambegaokar and Griffin ||5^ time or- 
dering operator Tr and time derivative operators dr do 
not commute and neglecting this subtlety can lead to for- 
mally divergent frequency sum matio ns [ p4|. Ta king heed 
of this subtlety, substitution of ( |C16 ) into(C14) amounts 
to 



dre 



4'l(l)vl/0(2)vi/t (3)vI/,(4) - ■^l{l)-^(,{2)^\{-i)i>s{A) 
-vl/t(l).i,^(2)vi;t(3)vi,,(4) + vl/t(i)vi,^(2)vl/t(3)vi,,(4)^, 

(C17) 



where the notation follows Eq. (B5 ) and Eq. ( p. 91 ). Upon 
utilizing the identities (|AE| ) and ( A13 ) and performing 
the standard Matsubara summation we have 

V^,{in) = ^Y.^n\V^,\rn){m\V''\n){e,, + e,r^fSr,r,m) 

nni 

(C18) 



where 5'mn is given by Eq. (A7). Analytically continuing 
iQ + iO we obtain 



-Hl7-|-iO 



(/n-/m). (C19) 



Note that the only difference between the Kubo contri- 
bution to t he thermal response (C19) and the spin re- 
sponse (B7) is the value o f the coupling constant. In the 
case of thermal response ( C19| ), the coupling constant is 
(cn + em)/2 which is eigenstate dependent, while in the 
case of spin response (B7) it is ?i/2 and eigenstate inde- 
pendent. 



Using Eq. (C12) we find that the Kubo contribution 
to the thermal transport coefficient is given by 



E 



(Cn + tmf 



iO) 



-\Tnni-\Tmn ( r 
2 f 



This can be written as 
where 



frn). 

(C20) 
(C21) 



(C22) 
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and 



(C23) 



Simil arly, we can separate the "diathermal" contribution 
( |CT^ as 



(C24) 



where M^^J and M^J refer to the first and second term 



in ( pl3 ) respectively. Using the completeness relation 
X^m 1"^) ('^1 = 1 it is easy to show that 



(C25) 



Comparison of ( |C23|) and (|C2l ) yields M^i^J+xj^J = 0. 
Therefore the thermal response coefficient is given by 



(C26) 



Utilizing commutation relationships (3.6), M^i^ can be 
expressed in the form: 

(C27) 

where the integral extends over the entire real line. The 
rest of the section follows closely Smrcka and Sti^eda . 
We define the resolvents G^: 



= {r] ± iO - Ho)-^ 



(C28) 



and operators 



B{r]) = zTr (v^G+V,6{f] - Ho) - V^5{tj - Ho)V,G- 

(C29) 

To facilitate Sommerfeld expansion we note that 



response coefficients kI]){T), M^^l'iT), af,^{T) have 
generic form 



(C30) 



L{T) = J f{rj)l{T^)drj, 
and after integration by parts 

Here L(^) is defined as 



Kv)d'n- 



(C31) 



(C32) 



Note that L{T = 0) = L(^ = 0) and in particular the spin 
conductivity at T = satisfies af^^{T^O) = ct^^(^ = 0). 
Identities (C31) and ( p32 ) will enable us to express 
the coefficients at finite temperature through the coeffi- 
cient s at zero temperature. For example, it follows from 
Eq. ( |C2^ ) that 



r^Tr (S{rj - Ho)ix^'V'' - x^V^") 



and from Eqs.(g8|, ^2g) 

^ J —00 



(C33) 



(C34) 



Similarly, coefficient Kj^J from (C22) can be written as 



so that 



dvf{v)v^S{v 



j/'nrriT/'Tnn 



iv - em + *0)2 (77 



(C35) 



J — 00 



,(?7-e™ + »0)2 (r,-e„-i0)2 
Using definitions ( |C29| ), i^^^^O 

can be expressed as 



(C36) 



r]^A{r])dr]+ / 7jB{r])drj. (C37) 



1 dB 
2~d^ 



drj. 



After integration by parts one obtains 

—00 — OC 

(C38) 

As shown in Ref . ( |Q ) the last term in this expression is 
exactly compensated by M'^^^(^): T his b ecomes evident 
after noting that definitions in Eq. (C2£) imply 



1 dB{r]) 

2 dr] 



-Tr 



dSjy - Ho) 
drj 



(C39) 



Substituting the last identity into the Eq.( |C38| ) and in- 
tegrating the second term by parts we obtain 



KilHO^e I A{^)drj 



r]Ti- (^d{r] - Ho){x''V'' ~ x^V^')'^ dT] (C40) 
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The second term here is equal to —M^tJ from (C33). 
After the cancellation the result simply reads: 

^?.(0 = 4'i(0+MW(e)=e' r (C41) 



Or, using (|C3J) 



<^%iO- (C42) 



Finally, from Eq.(C31) we find 
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